Superdiffusive comb: Application to experimental observation of anomalous diffusion 
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A possible mechanism of superdiffusion of ultra-cold atoms in a one-dimensional polarization 
optical lattice, observed experimentally in [Phys. Rev. Lett. 108, 093002 (2012)], is suggested. 
The analysis is based on a consideration of anomalous diffusion in a fractal comb [Phys. Rev. E 
83, 052106 (2011)]. It is shown that the transport exponent is determined by the fractal geometry 
of the comb due to recoil distributions resulting in Levy flights of atoms. 
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Recently, ID heavy-tailed distributions were observed 
in experimental studies of anomalous diffusion of ultra- 
cold 87 Rb atoms in a one-dimensional optical lattice 
It was found that the initial ensemble of atoms spread 
superdiffusively, such that the "full widths at half the 
maximum" (FWHM) increases with time like t~ with 
diffusion exponent 1 < [i < 2. Another important ob- 
servation was the dependence of the transport exponent 
on the depth of the lattice potential [lj. The theoreti- 
cal explanation of this fact, presented within the stan- 
dard semiclassical treatment of Sisyphus cooling [2j, is 
based on a study of the microscopic characteristics of the 
atomic motion in optical lattices and recoil distributions 
resulting in macroscopic Levy flights in space, such that 
the Levy distribution of the flights g(l) depends on the 
lattice potential depth Q: g(l) ~ l~ x ~ v with v = 
where D = cEn/Utra P , while [/trap is the lattice potential 
depth scaled by recoil energy En, and c is a dimensionless 
parameter (adopting notation from Ref. Q). A relation 
between the diffusion exponent /i and the Levy distri- 
bution exponent v was established for different regimes 
of the atomic dynamics, which is described by Fokker- 
Planck dynamics in an asymptotically logarithmic po- 
tential Hi]. 

Here, a relation between the transport exponent and 
the lattice potential depth is established in the framework 
of a fractional comb model [f| • We analyze an experimen- 
tal observation of dependence of a diffusion exponent as a 
function of the lattice depth [/trap, presented in Fig. 3 of 
Ref. [l[ (see here Fig. 1), where the limiting cases corre- 
spond to normal diffusion and ballistic motion, and these 
can be easily modelled by turbulent diffusion on a comb. 
This anomalous diffusion is described by the 2D distri- 
bution function P = P(x,y,t), and a special behavior 
is that the displacement in the ^-direction is possible 
only along the structure axis (a;-axis at y — 0), while any 
spreading along the x direction inside the fingers (motion 
along the y direction) is not possible. The Fokker-Planck 
equation in some dimensionless variables reads ( see e.g., 
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FIG. 1: A sketch of experimental observation of a wave packet 
spreading in a ID optical lattice [lj . The diffusion (fractional) 
exponent, presented vs the lattice potential depth, describes 
the full width at half the maximum (FWHM) of the atom 

distribution: FWHM ~ t» . The dashed lines correspond to 
the limiting cases. 



which can be defined by its Fourier image (see e.g., [g) 

^ l m] = -\^m (2) 

and P[f{x)] = f(k) with q = 1 in Eq. (J). It should 
be admitted that the y axis is the auxiliary space, in- 
troduced for the trap modelling. In other words, it is 
introduced to model a non-Markovian process by means 
of Markovian description. The true distribution is the 
distribution function along the x axis, which is 



P(x,t) 



P(x,y,t)dy. 



(3) 



dtP = 6(#)d\ x \P + d£P, 



(1) 



where d\ x \ is a particular case of the Riesz derivative, 



To establish the connection between turbulent diffu- 
sion on the comb with anomalous diffusion of cold atoms, 
we, first, model the limiting cases, shown in Fig.l. An 
effective constant diffusion coefficient in the y direction 
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can be considered as a function of the lattice depth 
d = d(U trap ). Therefore, Eq. ^ describes the two 
limiting cases (dashed lines in Fig. 1), where we use 
that d = 0(f7 t rap) is the Heaviside function. One easily 
checks [7j the full width at half the maximum (FWHM) 
of the atom distribution. For [/trap = one has FWHM 
= W (x 2 ) = t that corresponds to ballistic motion, while 
for d = 1 one finds normal diffusion with FWHM = ti . 
Here we use the definition 



(x 2 (t)) 



x 2 P(x, t)dx . 



(4) 



For the "shallow" trap lattice potentials, we modify 
Eq. ((T|) in the following way. From the experimental re- 
alization we know that after photon emission due to recoil 
atoms "fly" on distances distributed by power law and, 
correspondingly the fingers, as the traps, are distributed 
by power law with the fractal dimension v, related to the 
Levy flights. Therefore, this experimental realization can 
be described by the fractal comb model, developed in Q . 
Following this consideration, one can consider this set of 
the traps as a fractal set Fp(x) with the fractal dimen- 
sion which is embedded in the ID space, < v < 1. 
Therefore, the effective diffusion coefficient becomes in- 
homogeneous d — ¥ dx(x), where x( x ) is a characteristic 
function of Fo{x) 1 such that x( x ) = 1 f° r x S Fa(x) and 
x( x ) = for x (f. Fp(x). Taking this into account, we 
modify Eq. ([1} in the form 



d t P = S(y)d 2 x P + d X (x)d 2 v P. 



(5) 



To arrive at the corresponding modification of Eq (fT). 
we apply the Fourier transform to Eq. ([5]) with respect 
to the x coordinate. To this end we use the auxiliary 
identity 

/x poo 
x(y)f(y)dy = -d x / x(y)f(y)dy 
-oo J X 

(6) 

with the boundary conditions P(x = ±oo) = 0. This 
integration with the characteristic function can be carried 
out by means of a convolution [g, |8|, [9| 



x(y)f(y)dy => _«,.£/(*) = 



(7) 

where is the Gamma function and we also use the 
convenient notations of fractional integration: _ 00 I%f(x) 
S E3) El ■ We also used here the following arguments for 
the characteristic function. Note that 

/oo />oo 
X{y)f{y)dy= ^ / f(y)S(y-Xj)dy, 
-OO „ /— c 1 — OO 



where 



Xj) = M'(x) ~ |a;| c 



is a local fractal density, such that J^ x dM(y) ~ \x\ a cor- 
responds to the fractal volume. Therefore, due to Theo- 
rem 3.1 in Ref. [§| we have Jq f(y)dM.{y) ~ j^jj Jq(x — 

//:•" ;/':.'/:•<///. 

Using this fractional integration one obtains from Eq. 
© the fractional derivative dJ -oqI^P Qe, y, t)] defined in 
the Riemann-Liouville form [y, [ToL . Performing the 
Fourier transform and taking the symmetrical form, one 
obtains the following change in Eq. (p} d — > d|/c| 1_iy . 
Note that v = £(C/trap) is the fractal dimension of the 
fingers distribution on the x axis. After this change Eq. 
([!]) reads 



d t P = -5{y)\k\P + d\k\ 2 ~ v d 2 y P , 



(8) 



where v = 1 + v corresponds to the Levy distribution 
related to the lattice potential depth Q. To satisfy the 
limiting cases, we have v = 2 and d = Q(U trap ) is the 
Heaviside function. 

Our aim now is to find P(x,t), defined in Eq. ([3]) 
for "shallow" trap lattice potentials, when 1 < v < 2 
Q < D < |). To this end, an analysis developed in Q 
is applied. One carries out the Laplace transform in the 

time domain £[P(k, y,t)} = P(k,y,s). Looking for the 
solution of the Laplace image in the form 



P(k, y, s) - exp[-|y|v/|fc|- 2 S /d]/(fc, a) . 



(9) 



one arrives at the intermediate expression in the form of 
the Laplace and Fourier inversions 



P{x,y,t)=F- 1 lt 



2e 



\vWs\k\—*/d 



2^Jsd\k\ 2 ~» + \k\ 



(10) 



Integration over y and the inverse Laplace transform 
yield a solution in the Fourier inversion form: 



P(x,t) 
Here 
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-^£ h [- l -^WIJd)dk. (11) 



£ a (-z) 



e u du 



2m ./ 7 u a + z 



is the Mittag-Leffler function defined by the inverse 
Laplace transform with a corresponding deformation of 
the contour of the integration [12j . First, we admit the 
scaling variable x/t u that corresponds to the superdiffu- 
sion expansion 

FWHM - ti . 

This also corresponds to the experimental observations, 
presented in Ref. [l|, and to the scaling obtained in 
Ref. (see Eq. (11) there). 

Let us consider an asymptotic behavior of this su- 
perdiffusion expansion of the initial packet of ultra-cold 
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atoms at x — > oo (k —> 0). In this case, the argument of 
the Mittag-Leffler function is small, yielding [10] 

fa(z) ^ ex K-f(TT^))' 

and one obtains 

^■ t )-s£«-*-»(-^)*- 

This is nothing but the Fourier inversion of the the char- 
acteristic function of a centered and symmetric Levy dis- 
tribution [ljj, E3| that describes Levy flights of atoms. 
This integration yields the analytical solution in the 
form of the Fox function with the power-law asymptotics 

MM 

? (s>*)~£jT+^ v < 2 - ( 12 ) 

From here one obtains for the first moment [lol [13} 
(x(t)} ~ t~ , which corresponds to the scaling obtained 
above. 



In conclusion, we note that this description is a partic- 
ular case of a general scheme of the wave-packet spread- 
ing described by the fractional Fokker-Planck equation 
(FFPE) that can be obtained for the true distribution 
(fT2j). Integrating the Laplace image of Eq. (|5|) over 
y, and taking into account Eq. (j9]), one performs the 
Laplace and the Fourier inversions that yields the FFPE 
for the true distribution 

d}p[x,t) = ^j=dlf{x,t). (13) 

Here 9 t 2 is the Caputo time fractional derivative [T3| and 
d£l is the Riesz fractional derivative [6j] defined in Eq. 
This equation is a particular case of the FFPE 

d?P = K(a,d)d^P, (14) 

where K(a,d) is a generalized diffusion coefficient and 
the scaling q/a = v is fulfilled for the initial ensemble 
of atoms, which spreads superdiffusively with a < 1 and 
q < 2. 
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